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(54) Channel allocation in radio systems. 



(57) i A method of allocating radio channels, /be- 
— tj^errth^ase^tatioris of -a ■m u l t i -base s ta tion 
radio system, uses a mathematical model of the 
radio system which is derived from an analog of 
the radio system. The analog is a thermal bath 
containing a plurality of interacting particles, 
each particle in a particular spin state. The 
number of particles at a given location in the 
thermal bath, corresponding to the traffic de- 
mand at an equivalent location in the radio 
system, each occupied spin state correspond- 
ing to a particular radio channel. The model is 
operated by allowing the spin states to relax, by 
means of a stochastic reduction in temperature 
of the thermal bath, from an initial infinite tem- 
perature, to a minimum energy, which equates 
to an acceptable level of interference between 
radio channels, operated in the different cells of 
the radio system. The distribution of spin states 
at the minimum energy determines the allo- 
cation of radio channels. Apparatus for imple- 
menting the method includes a central control 
unit which runs the mathematical model and 
issues commands to the individual base sta- 
tions allocating channels. 
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The invention relates to a method of allocating ra- 
dio channels between the base stations of a radio 
system, having a plurality of base stations, and appa- 
ratus for implementing the method. 

In radio systems used for telecommunications, 
especially, but not exclusively, those available for use 
by the general public, such as mobile radio telephony 
systems, traffic demand varies according to the time 
of day, from day to day, and with the season of the 
year. This phenomena is a problem with which all 
those engaged in the arts of operating, installing and 
designing such systems are familiar. In the first multi- 
base station radio systems, channel allocation was 
implemented when the radio systems were initially 
commissioned. Prior to commissioning measure- 
ments relating to field strengths, ranges of coverage, 
interference characteristics etc. were carried out. Re- 
distribution of radio channels was carried out when 
new base stations and(or new channels were intro- 
duced to the systems. Redistribution of channel allo- 
cation was also carried out if it was found that the traf- 
fic capacity of the systems was inadequate. The util- 
ity of such radio systems remained constant for long 
periods because changes to channel allocation re- 
quired a significant amount of work to be performed 
on the systems. 

In multi-base station systems of the type descri- 
bed above it was discovered, as traffic levels in- 
creased, that traffic demand was concentrated on 
certain geographic areas. Furthermore these geo- 
graphic areas changed with time, so that it became 
desirable to redistribute available radio channels be- 
tween base stations. 

Initially, channel redistribution was carried out by 
mapping the traffic demand at different times of day. 
This information lead to the conclusion that channel 
distribution should be changed at certain times of the 
day. However, it was soon found that this solution was 
too clumsy to be functional. To overcome this problem 
those engaged in the field of radio communications 
attempted to set up, in a variety of ways, automatic 
systems for distributing radio channels, between 
base stations, in an optimal manner. 

For a proper analysis of traffic demand knowl- 
edge of the following factors is required: 

- the geographic location of base stations; 

- frequencies of the available channels; 

- the transmitter powers utilised; 

- interactions of the channels with each other in 
different circumstances; etc. 

A variety of methods have been used for deter- 
mining the optimal redistribution of channels between 
base stations, for example linear programming with a 
number of subsidiary conditions, reflecting the vari- 
ables referred to above. However, one difficulty is 
that the formulation of the problem is particularly 
complex, and requires a computer aided solution. 
Computers are indispensable to today's flexible radio 



systems. 

The prior art methods h\\ have the disadvantage 
that the mathematical mooeis which they employ be- 
come unmanageable with increasing numbers of 

5 parameters and conditions. The solution to this prob- 
lem is of comparable difficulty to that of the well 
known "travelling salesman" problem. 

Furthermore it is not feasible, in practice, to carry 
out field measurements every time the utilisation of 

10 radio channels is to be changed. 

The present invention relates to a method of de- 
termining the channel allocation for multi-base station 
radio systems which can be used in cases where, be- 
cause of the complexity of the problem, the earlier 

is me. hods cannot be employed. Apparatus for imple- 
me-uing the method is also described. A radio net- 
work which incorporates the present invention in- 
cludes means for collecting information needed to 
make decisions on channel allocation, together with 

20 a central control device which calculates the channel 
allocation, as between different base stations, and al- 
locates channels accordingly. The calculation of the 
channel allocation is made by reference to an analo- 
gous physical system which behaves in an analogous 

25 manner to a multi-base station radio system, namely 
a heat bath containing interacting particles. One 
group of particles represents a given allocation of 
channels to a given base station. By minimising the 
total energy of the system, the optimal channel allo- 

30 cation between base stations is determined. 

According to a first aspect of the present inven- 
tion there is provided a method of distributing chan- 
nels between base stations in a multi-base station ra- 
dio system by constructing a mathematical model 

35 from which can be determined a channel distribution 
which has an acceptable level of inter-channel inter- 
ference, and assigning channels to the base stations 
in accordance with the predictions of the mathemat- 
ical model characterised in that the mathematical 

40 model is based on a heat bath having a plurality of 
particles immersed therein, in which 

- each base station of the radio system corre- 
sponds to one particle, 

- each radio channel corresponds to a discrete 
45 state of said particle, 

- selected channels are represented by different 
states of a particle so that changes in said dis- 
tribution of channels are represented by 
changes in the states of the particles, 

so - interference between base station channels 

are represented by the interaction between the 
states of different particles, 

- a change in the interference range coverage 
area ratio (C/l) is represented by a change in 

55 the energy of the thermal bath, 

- the number of particles and states in the math- 
ematical model represents the total number of 
base stations and channels in the radio system 



BNSDOCID: <EP 0565499A1 J_> 




EP 0 565 



which corresponds to that necessary to meet 
a predetermined traffic demand with the max- 
imum acceptable level of interference between 
channels, 5 
a first condition of the radio system is specified by the 
number of particles and the interaction between them 
in the heat bath, a redistribution of the radio channels 
or a change in the number of channels is represented 
in the model by a change in the states and interac- 10 
tions of the particles, the channel allocation is estab- 
lished by setting an initial temperature for the heat 
bath at infinity and then reducing the temperature in 
steps so that the states of the particles change in ac- 
cordance with a stochastic process, and the energy 15 
of the heat bath is changed in order to achieve an en- 
ergy content which corresponds to an acceptable in- 
teraction between particles and acceptable level of 
interference between base stations. 

According to another aspect of the present inven- 20 
tion there is provided, for a multi-base station radio 
system having a plurality of base stations, each base 
station having one or more radio channels allocated 
to it, a method of distributing available radio channels 
between base stations so as to reduce interference 25 
between radio channels allocated to different base 
stations to an acceptable level is characterised in that 
a distribution of radio channels between base stations 
is determined from a mathematical model based on 
a plurality of particles in a thermal bath, in which par- 30 
ticle(s) at a given location in the thermal bath corre- 
spond to a given base station, and the particle(s) 
have a plurality of states, each occupied state corre- 
sponding to a radio channel, the occupancy of states 
determined by a thermal distribution, and the partic- 35 
les having interaction energies corresponding to in- 
terference between base stations, a desired alloca- 
tion of radio channels being determined by setting the 
initial temperature of the thermal bath at infinity and 
then reducing the temperature of the thermal bath in 40 
steps so that a local minimum energy is attained by a 
stochastic process in which the interaction between 
particles corresponds to an acceptable level of inter- 
ference between base stations and determining the 
radio channel allocation from the distribution of occu- 45 
pied states of the particle(s). 

According to a further aspect of the present in- 
vention there is provided in a multi-base station radio 
system having a plurality of base stations, each base 
station having one or more radio channels allocated so 
to it, apparatus for distributing available radio chan- 
nels between base stations so as to reduce interfer- 
ence between radio channels allocated to different 
base stations to an acceptable level, includes a con- 
trol means for sending commands to individual base 55 
stations determining the radio channels allocated to 
the base station, which control means is arranged so 
that it can receive data on traffic demand in the radio 
system and the availability and utilisation of radio 
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channels, is characterised in that the control means 
includes a computer on which is run a mathematical 
model of the radio system from which a distribution of 
radio channels between base stations is determined, 
the mathematical model based on a plurality of par- 
ticles in a thermal bath, in which particle(s) at a given 
location in the thermal bath correspond to a given 
base station, and the particle(s) have a plurality of 
states, each occupied state corresponding to a radio 
channel, the occupancy of states determined by a 
thermal distribution, and the particles having interac- 
tion energies corresponding to interference between 
base stations, a desired allocation of radio channels 
being determined by setting the initial temperature of 
the thermal bath at infinity and then reducing the 
temperature of the thermal bath in steps so that a lo- 
cal minimum energy is attained by a stochastic proc- 
ess in which the interaction between particles corre- 
sponds to an acceptable level of interference be- 
tween base stations and determining the radio chan- 
nel allocation from the distribution of occupied states 
of the particle(s). 

According to yet another aspect of the present in- 
vention there is provided a multi-base station radio 
system having a plurality of base stations, each base 
station having one or more radio channels allocated 
to it, including apparatus for distributing available ra- 
dio channels between base stations so as to reduce 
interference between radio channels allocated to dif- 
ferent base stations to an acceptable"! evel, which ap- 
paratus includes a control means for sending com- 
mands to individual base stations determining the ra- 
dio channels allocated to the base station which con- 
trol means is arranged so that it can receive data on 
traffic demand in the radio system and the availability 
and utilisation of radio channels, is characterised in 
that the control means includes a computer on which 
is run a mathematical model of the radio system from 
which a distribution of radio channels between base 
stations is determined, the mathematical model 
based on a plurality of particles in a thermal bath, in 
which particle(s) at a given location in the thermal 
bath correspond to a given base station, and the par- 
ticle^) have a plurality of states, each occupied state 
corresponding to a radio channel, the occupancy of 
states determined by a thermal distribution, and the 
particles having interaction energies corresponding 
to interference between base stations, a desired al- 
location of radio channels being determined by set- 
ting the initial temperature of the thermal bath at in- 
finity and then reducing the temperature of the ther- 
mal bath in steps so that a local minimum energy is 
attained by a stochastic process in which the interac- 
tion between particles corresponds to an acceptable 
level of interference between base stations and de- 
termining the radio channel allocation from the dis- 
tribution of occupied states of the particle(s). 

Embodiments of the present invention will now be 
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described by way of example with reference to the ac- 
companying drawings in which: 

Figure 1 illustrates diagrammatically a mobile ra- 
dio network with a plurality of base stations whose in- 
terference ranges completely overlap, together with 
an analogous heat bath in which particles are im- 
mersed. 

Figure 2 shows a mobile radio network of the 
same type as that shown in Figure 1, in which the in- 
terference ranges have been reduced to a number of 
islands, together with an analogous particle contain- 
ing heat bath. 

Figure 3 shows a mobile radio network of the 
type shown in Figure 1 , in which the interference level 
is acceptable, together with the analogous particle 
containing heat bath. 

Figure 4 shows a base station with channels k1 , 

k2 kn, together with particle(s) which represents 

this base station in the associated heat bath. 

Figure 5 shows a schematic representation of a 
radio system according to the invention. 

A radio system 1, (see Figure 5), consists of a 
number of radio base stations 2, 3, 4, which commu- 
nicate with mobile transceivers 8, 9, 11, within cells 
32, 33, 34. Each base station 2, 3. 4, is allocated one 
or more radio channels. The traffic demand varies at 
different times of the day, and when this occurs there 
is a need to vary the channel utilisation. It may also 
be necessary to rearrange the channel allocation 
when new stations are brought into the system, or sta- 
tions are taken out of the system. 

When traffic demand changes, a command 
which causes a rearrangement of the channel alloca- 
tion is initiated from a central control unit 20, in the 
system. The central control unit includes a computer 
21. The traffic demand via each base station 2, 3, 4, 
is determined together with the restrictions relating to 
permissable channel selections. Examples of restric- 
tions are: 

- the channel demand via each base station; 

- the number of possible channels which can be 
handled by each base station; 

- which channels can be adjacent to each other 
in the same base station; and 

- which channels can be adjacent to channels al- 
located to neighbouring base stations. 

Because a thermal bath 10, containing particles 
is analogous to a multi-base station radio system 1, 
each base station 2, 3, 4, with its associated channels 
and restrictions can be treated theoretically as a 
group of particles 13,14,15, each in a given quantum 
spin state, in a thermal system. Each particle inter- 
acts with the other particles in the system. The quan- 
tum spin states of a group of particles, in the thermal 
system, corresponds to a given base station with its 
allocated radio channels, illustrated diagrammatically 
in Figure 4. A given particle spin state in the thermal 
bath analog corresponds to a given channel in the ra- 
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dio system. The traffic demand in a given cell of the 
radio system corresponds to the number of particles 
at the corresponding location in the thermal bath 

5 analog. The different restrictions on channel alloca- 
tion can be given a weighting, depending on their sig- 
nificance, for controlling the system's traffic capacity 
and usability. Such restrictions can be given effect in 
the model through exclusion principles on occupancy 

10 of certain states and probabilities of occupancy of 
certain states. 

Variations on this model can easily be envisaged, 
for example, vibration states of particles could be 
used instead of spin states of particles. Furthermore 

15 a base station could be represented by a single par- 
ticle having a plurality of spin states a number of 
which can be simultaneously occupied c.f. an atom 
with multiple electronic spin states. The key element 
of the model is that the individual particle(s) repre- 

20 senting base stations must have a series of discrete 
states which can be thermally excited and are capa- 
ble of interacting with the states in other particles. 

A variety of techniques are available for mathe- 
matically modelling a collection of interacting partic- 

25 les, having different spin states, in a thermal bath and 
the reader is referred to the following text books for 
details of the available techniques: 

- C. Kittel. 'Thermal Physics", pub. John Wiley & 
Sons Inc. 1969; 

30 - P.G. Moel. S.C. Port and C.J. Stone, "Introduc- 

tion to Stochastic Processes", pub. Moughton 
Mifflin Company, 1972; 

- S. Brunak and B. Lantrup, "Neural Networks", 
pub. World Scientific, 1990; and 

35 - E. Aarts and J. Korst, "Simulated Annealing 

and Boltzmann Machines", pub. John Wiley & 
Sons Inc., 1989 
Whatever technique is used to model the analo- 
gue system, the thermal model is initially given an in- 

40 finite, or very high, initial temperature, see Figure 1. 
In this state the interference ranges of the corre- 
sponding mobile radio system overlap completely, 
and in the analogue thermal bath the particles inter- 
act strongly. Thereafter, lower temperatures are simu- 

45 lated. Initially islands of interference 5, 6, 7, are 
formed, see Figure 2. Eventually the temperature is 
reduced to the point at which the islands of interfer- 
ence 5, 6, 7, have contracted to an acceptable level, 
see Figure 3. In this calculation, interference ranges 

so in the radio system correspond to energies. As this 
process of progressively lowering the temperature is 
continued, a minimisation of the total system energy 
is sought. The minimum energy thus found corre- 
sponds to the channel distribution which produces 

55 the best utilisation of the system for a given traffic de- 
mand. The object is to arrive at the lowest local sys- 
tem energy which permits a given, minimum, degree 
of particle interaction, by allowing the different distrib- 
utions of states for the different particles to vary with 

4 
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temperature. 

In the model, the restrictions are determined from 
measured or model-calculated values as the ratio be- 
tween connection level and interference level (C/l) for 5 
each cell, cell pair and channel spacing. Interferences 
are obtained both from base stations 2, 3, 4, and mo- 
biles 8, 9, 10, in the reception areas (cells) 32, 33, 34. 
In many cases, it is adequate to consider co-channel 
and adjacent channel interferences. By co-channel 10 
interference is meant interference between the same 
channel used in different cells. Adjacent channel in- 
terference is interference between neighbouring 
channels operated in different cells. If the (C/l) ratio 
is lower than a minimum acceptable value, this ex- 15 
eludes the use of channels, within the relevant chan- 
nel spacing, from this cell, if the channel is utilised in 
a second adjacent cell. The limit value is determined 
by the probability that interference will occur. In some 
cases it is desirable to minimise the total interference 20 
level. In such cases exclusions are introduced which 
correspond to a higher value of C/l, or lower interfer- 
ence probability, but with a lower weight so that all 
necessary exclusions must still be fulfilled. In this 
case the value of the weights must be calculated by 25 
means of testing. To directly translate C/l and interfer- 
ence probabilities into weights, the optimisation func- 
tion (the energy) which makes up the sum of all 
weights corresponding to exclusions which have not 
been fulfilled, must be changed. Exclusions can also 30 
be obtained with respect to the system as a whole. 
The need for field measurements can be significantly 
reduced by estimating the co-channel exclusions 
from adjacent channel exclusions. 

The problem thus formulated can be thought of 35 
as a set of thermal equations with associated restric- 
tions which must be minimised. In the calculations, it 
is assumed that the thermal process follows a sto- 
chastic process, for example a Markov process. New 
conditions in the system are simulated by dropping 40 
the temperature in the analogue thermal system, see 
Figures 1 to 3. In particular Figure 2 shows a radio 
network 1, in which the interference ranges are re- 
duced to a number of islands 5, 6, 7. The interference 
level is still not acceptable and so in the thermal ana- 45 
logue the temperature must be lowered still further. 
Figure 3 shows a radio network in which the interfer- 
ence levels are acceptable. In the thermal analogue 
the temperature has been lowered to the point at 
which an acceptable interaction exists between the 50 
particles. This implies that the radio system as a 
whole is functioning in an acceptable manner. The 
temperature is dropped by lowering the energy con- 
tent in the system. It is important in this context that 
sub-optimisation is avoided, since this involves local 55 
minima being utilised in the final solution without ach- 
ieving the best solution for the system. This can be 
avoided by reducing the temperature in energy steps 
which correspond to approximately one standard de- 



viation, of the energy distribution between particles 
present in the thermal bath. Selected temperatures 
are extrapolated from a curve which is adapted to the 
mean energies. 

The established restrictions may be exceeded 
during the course of the calculations, which leads to 
infringement of the stipulated rules for the operation 
of the radio system. Every time a change of energy 
is introduced, there either is or is not an infringement 
of the established restrictions. If there is no infringe- 
ment, the result is accepted. When an infringement of 
the established rules is observed, this implies that an 
energy addition is obtained in the system, which in 
turn implies that the risk of interference between 
base stations and(or radio channels is increased. De- 
pending on the importance of the deviation, the dif- 
ferent infringements are weighted which enables a 
better control over the end result. The weightings 
used in the system can also be changed, for example 
at different times of the day. 

A number of iterations are made up to an empiri- 
cal factor multiplied by a nominal relaxation time 
which comprises the ratio between the variance of 
the energy and the variance of the energy changes. 
To determine the variance of the energy changes, ex- 
trapolation is used. This does not prevent the use of 
other known methods for calculating the value. 

A temperature increase may sometimes be re- 
quired in difficult cases when the energy is near op- 
timum. Under such circumstances, it may even be 
necessary, occasionally, to concentrate the calcula- 
tions on small, but especially troublesome, subre- 
gions, where, normally, standard methods suffice for 
channel allocation. It is also possible to allow the tem- 
perature to vary somewhat over the system depend- 
ing on. among other things, weightings. 

The following points should be noted about the 
thermal bath analog and its relationship to the radio 
system which it simulates: 

- The particles in the thermal bath are arranged 
in groups at locations corresponding to the 
base stations in the radio system; 

- Each particle has a discrete or quantum spin 
state corresponding to a given channel in the 
radio system; 

- The number of states corresponds to the num- 
ber of channels, however the energy of a par- 
ticular state at a particular location is deter- 
mined by the interaction with other spins (i.e. 
interference) rather than by any external field. 
However an external field, such as a magnetic 
field could be introduced if, for some reason, 
certain channels are preferred to others; 

- The particles in the thermal bath are in "quasi- 
equilibrium" with each other, that is to say they 
are in equilibrium in a state space which con- 
tracts as the bath is cooled, excluded states 
simply become too improbable to occur; 
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- Interference between channels corresponds to 
interactions between the spin states, and con- 
tributes to the overall heat capacity of the bath; 

- The traffic demand in a cell corresponds to the 
desired number of particles at the correspond- 
ing location in the thermal bath. Usually the 
number of particles at a given location is kept 
constant during each simulation, although 
there is no overriding requirement for the num- 
ber of particles at each location to be kept con- 
stant; 

- The transmitter power of a given base station 
is one of the parameters which determine the 
strength of interaction between locations; 

- The C/l ratio corresponds to the heat capacity 
of the system and relates to the system energy 
through the heat capacity, the temperature is 
the controlling parameter and is frequently kept 
uniform across the whole system; 

- As more channels are added to a base station, 
this can be regarded as the creation of more 
states, or as an increased occupancy of the ex- 
isting states. 

System constraints on the radio system can be 
given effect in the model through the use of exclusion 
principles, for example preventing multiple occupancy 
of the same state at a given location (equivalent to a 
base station being unable to make multiple use of the 
same channel), and by means of prohibited states at 
certain locations. 

In normal use the number of channels allocated 
to each base station is fixed, the model determines 
not the number of channels allocated to a given base 
station, but which channels are allocated to a given 
base station. It is however possible to vary the num- 
ber of channels allocated to a base station if this is 
really necessary, this corresponds to changing the 
number of particles at a particular location in the ther- 
mal bath. 

As the temperature of the thermal bath is re- 
duced the spin states of the individual particles will re- 
arrange themselves to produce a lower energy con- 
figuration. This takes into account the interaction en- 
ergies between the particles, which make a contribu- 
tion to the heat capacity of the system. Eventually a 
temperature will be reached which corresponds to a 
local minimum which, hopefully, gives acceptable in- 
terference levels. Because of the presence of multi- 
ple local energy minima, the final minimum achieved 
is not necessarily the global minimum. If the interfer- 
ence level is not satisfactory it may be necessary to 
raise the temperature of the thermal bath in order to 
hunt for a new local minimum. 

The final allocation of channels is read from the 
model by determining the spin states of the particles 
at each location. 
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EXAMPLE 

It is by no means easy to visualise just how the 
5 analogy between a multi-base station radio system 
and a thermal bath containing particles operates in 
practice. A rather simple and trivial example can be 
constructed, and is presented here, not as an exam- 
ple of an actual system, but as an illustration of how 
10 the mathematical model used in the invention can be 
constructed. 

Consider a radio system with three base stations 
A, B, C and two available channels 1. 2. 

From field strengths and(or model computations. 
15 the interaction strengths for adjacent channel inter- 
ference are determined to be: 
A - B : 0.4 
B - C : 0.3 
A-C:0.1 

20 From field strengths and(or model computations, 

the interaction strengths for co-channel interference 
are determined to be: 
A- B : 1.0 
B - C : 1.0 
25 A - C : 0.5 

The traffic demand i.e. the number of channels 

is: 

A : 2 
B : 1 

30 C : 1 

Site A is represented by a vector V a = [a 1t a 2 ], 
where a ; = 0 or ai = 1 , a^ + a 2 = 2 (the traffic demand). 
Sites B and C are similarly represented. 
The energy of the system is thus 
35 E = {1.0 * (a, * b, + a 2 * b 2 ) + 1.0 * (b, * + 
b 2 * c 2 ) + 0.5 * (a! * c, + a 2 * c 2 ) + 0.4 * (a^ * b 2 
+ a 2 * bi) + 0.3 * (bi * c 2 + .b 2 * c,) + 
0.1 * (a, * b 2 + a 2 * bO} 
The vectors V a , V b , and V c are initialised in a ran- 
40 dom fashion, i. e. V a = [1,1], V b = [0,1], V c = [1,0]. 
(There is only one possibility for V a in this example, 
but we will ignore this). 

The state transitions are generated via a Markov 
chain, so that the probability for a state E is given by 
45 a thermal distribution, which is usually the Boltzmann 
distribution P a exp (-E/T) where T is the temperature. 
The temperature is high (actually infinite) at the be- 
ginning and then gradually lowered. Finally the sys- 
tem will settle into a local energy minimum which, in 
50 general, is close to a global minimum, thus minimizing 
interference. 

The method of the present invention has the ad- 
vantage over prior art methods of enabling much 
more complex forms of the problem of channel allo- 
55 cation to be solved. 

The results obtained from the calculations are 
used for the automatic redistribution of channels in a 
radio system, which redistribution is effected by a 
central control unit, which allocates the available 
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channels in the best manner at a given time. 



Claims 5 

1. A method of distributing channels between base 
stations (2, 3, 4) in a multi-base station radio sys- 
tem (1) by constructing a mathematical model 
from which can be determined a channel distrib- 10 
ution which has an acceptable level of inter-chan- 
nel interference, and assigning channels to the 
base stations in accordance with the predictions 
of the mathematical model characterised in that 
the mathematical model is based on a heat bath 15 
(10) having a plurality of particles (13, 14, 15) im- 
mersed therein, in which 

- each base station (2, 3, 4) of the radio sys- 
tem (1) corresponds to one particle, 

- each radio channel corresponds to a dis- 20 
crete state of said particle, 

- selected channels are represented by dif- 
ferent states of a particle so that changes in 
said distribution of channels are represent- 
ed by changes in the states of the particles 25 
(13, 14, 15), 

- interference between base station chan- 
nels are represented by the interaction be- 
tween the states of different particles, 

- a change in the interference range cover- 30 
age area ratio (C/l) is represented by a 
change in the energy of the thermal bath 
(10), 

- the number of particles and states in the 
mathematical model represents the total 35 
number of base stations (2, 3, 4) and chan- 
nels in the radio system (1) which corre- 
sponds to that necessary to meet a prede- 
termined traffic demand with the maximum 
acceptable level of interference between 40 
channels, 

a first condition of the radio system (1) is speci- 
fied by the number of particles (13, 14, 15) and 
the interaction between them in the heat bath 
(10), a redistribution of the radio channels or a 45 
change in the number of channels is represented 
in the model by a change in the states and inter- 
actions of the particles (13. 14, 15), the channel 
allocation is established by setting an initial tem- 
perature for the heat bath (1 0) at infinity and then so 
reducing the temperature in steps so that the 
states of the particles (13, 14, 15) change in ac- 
cordance with a stochastic process, and the en- 
ergy of the heat bath (10) is changed in order to 
achieve an energy content which corresponds to 55 
an acceptable interaction between particles (13, 
14, 15) and acceptable level of interference be- 
tween base stations (2, 3, 4). 



2. A method as claimed in claim 1 further character- 
ised in that the temperature and energy of the 
heat bath (10) is allowed to drop in steps selected 
with respect to the mean temperature. 

3. A met hod as claimed in claim 2 further character- 
ised in that the temperature of the heat bath (10) 
is reduced in energy steps which correspond to 
one standard deviation of the energy distribution 
from the mean energy in order to avoid sub-op- 
timisation. 

4. A method according to either claim 2 or 3 further 
characterised in that the temperature drops are 
selected with respect to the mean energy accord- 
ing to a prediction method. 

5. A method according to any preceding claim fur- 
ther characterised in that a number of iterations 
is automatically selected so that improved con- 
vergence is obtained at low energies. 

6. In a multi-base station radio system (1) having a 
plurality of base stations (2, 3, 4) each of which 
has one or more radio channels allocated to it, ap- 
paratus for allocating radio channels to the base 
stations having control means (20) for effecting a 
redistribution of channels between base stations 
(2, 3, 4) which control means (20) receives data 
on traffic demand in the radio system (1) and the 
availability and utilisation of radio channels, is 
characterised in that the control means (20) in- 
corporates a mathematical model for the deter- 
mination of channel distribution between base 
stations (2, 3, 4) which operates according to the 
method claimed in any one of claims 1 to 5. 

7. A multi-base station radio system (1) having a 
plurality of base stations (2, 3, 4) each of which 
has one or more radio channels allocated to it, in- 
cluding apparatus for allocating radio channels to 
the base stations (2, 3, 4) having control means 
(20) for effecting a redistribution of channels be- 
tween base stations (2, 3, 4) which control means 
(20) is adapted to receive data on traffic demand 
in the radio system (1) and the availability and 
utilisation of radio channels, is characterised in 
that the control means (20) incorporates a math- 
ematical model for the determination of channel 
distribution between base stations which oper- 
ates according to the method claimed in any one 
of claims 1 to 5. 

8. For a multi-base station radio system (1) having 
a plurality of base stations (2, 3, 4), each base 
station having one or more radio channels allo- 
cated to it, a method of distributing available radio 
channels between base stations (2, 3, 4) so as to 
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reduce interference between radio channels allo- 
cated to different base stations (2, 3, 4) to an ac- 
ceptable level is characterised in that a distribu- 
tion of radio channels between base stations (2, 5 
3, 4) is determined from a mathematical model 
based on a plurality of particles (13, 14, 15) in a 
thermal bath (10), in which particle(s) at a given 
location in the thermal bath (10) correspond to a 
given base station, and the particle(s) (13,14,15) 1 o 
have a plurality of states, each occupied state 
corresponding to a radio channel, the occupancy 
of states determined by a thermal distribution, 
and the particles having interaction energies cor- 
responding to interference between base sta- is 
tions (2, 3, 4), a desired allocation of radio chan- 
nels being determined by setting the initial tem- 
perature of the thermal bath (10) at infinity and 
then reducing the temperature of the thermal 
bath (1 0) in steps so that a local minimum energy 20 
is attained by a stochastic process in which the in- 
teraction between particles (13, 14, 15) corre- 
sponds to an acceptable level of interference be- 
tween base stations (2, 3, 4) and determining the 
radio channel allocation from the distribution of 25 
occupied states of the particle(s). 

9. A method as claimed in claim 8 further character- 
ised in that 

-. the states of the particles (13, 14, 15) are 30 
spin states, 

- each base station (2, 3, 4) is represented by 
a group of particles (13, 14, 15), 

- the number of particles in a group corre- 
sponds to the number of channels assigned 35 
to a base station represented by that group 

of particles, and 

- the spin state of a given particle in a group 
represents a particular radio channel. 

40 

10. A method as claimed in claim 8 or 9 further char- 
acterised in that the strength of interaction be- 
tween states of the particles (13, 14, 15) is deter- 
mined from field measurements on the radio sys- 
tem (1) and takes into account such factors as 45 
transmitter strength and local geography. 

11. A method as claimed in claim 10 further charac- 
terised in that the temperature of the thermal 
bath (10) is reduced in steps of approximately so 
one standard deviation of the mean thermal en- 
ergy of the particles (13, 14, 15) in the thermal 
bath (10), until an energy minimum is attained. 

12. A method as claimed in any one of claims 8 to 11 55 
further characterised in that the stochastic proc- 
ess is a Markov process. 

13. A method as claimed in any one of claims 8 to 12 
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further characterised in that the thermal distrib- 
ution is a Boltzmann distribution. 

14. In a multi-base station radio system (1) having a 
plurality of base stations (2, 3, 4), each base sta- 
tion having one or more radio channels allocated 
to it, apparatus for distributing available radio 
channels between base stations (2, 3, 4) so as to 
reduce interference between radio channels allo- 
cated to different base stations (2, 3, 4) to an ac- 
ceptable level, includes a control means (20) for 
sending commands to individual base stations (2, 
3, 4) determining the radio channels allocated to 
the base station (2, 3, 4), which control means 
(20) is arranged so that it can receive data on 
traffic demand in the radio system (1) and the 
availability and utilisation of radio channels, is 
characterised in that the control means (20) in- 
cludes a computer (21) on which is run a mathe- 
matical model of the radio system from which a 
distribution of radio channels between base sta- 
tions (2, 3, 4) is determined, the mathematical 
model based on a plurality of particles in a ther- 
mal bath (10), in which particle(s) (13, 14, 15) at 
a given location in the thermal bath (10) corre- 
spond to a given base station, and the particle(s) 
(13, 14, 15) have a plurality of states, each occu- 
pied state corresponding to a radio channel, the 
occupancy of states determined by a thermal dis- 
tribution, and the particles (13, 14, 15) having in- 
teraction energies corresponding to interference 
between base stations (2, 3, 4), a desired alloca- 
tion of radio channels being determined by set- 
ting the initial temperature of the thermal bath 
(10) at infinity and then reducing the temperature 
of the thermal bath (10) in steps so that a local 
minimum energy is attained by a stochastic proc- 
ess in which the interaction between particles 
corresponds to an acceptable level of interfer- 
ence between base stations (2, 3, 4) and deter- 
mining the radio channel allocation from the dis- 
tribution of occupied states of the particle(s) (13, 
14, 15). 

15. Apparatus as claimed in claim 14 further charac- 
terised in that 

- the states of the particles (13, 14, 15) are 
spin states, 

- each base station (2, 3, 4) is represented by 
a group of particles (13, 14, 15), 

- the number of particles in a group corre- 
sponding to the number of channels as- 
signed to the base station represented by 
that group of particles, and 

- the spin state of a given particle in a group 
represents a particular radio channel. 

16. A multi-base station radio system (1) having a 
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plurality of base stations (2, 3. 4), each base sta- 
tion having one or more radio channels allocated 
to it, including apparatus for distributing available 
radio channels between base stations (2, 3. 4) so 5 
as to reduce interference between radio chan- 
nels allocated to different base stations (2, 3, 4) 
to an acceptable level, which apparatus includes 
a control means (20) for sending commands to in- 
dividual base stations (2, 3, 4) determining the ra- 10 
dio channels allocated to the base station which 
control means (20) is arranged so that it can re- 
ceive data on traffic demand in the radio system 
(1) and the availability and utilisation of radio 
channels, is characterised in that the control 15 
means (20) includes a computer (21) on which is 
run a mathematical model of the radio system (1 ) 
from which a distribution of radio channels be- 
tween base stations (2, 3, 4) is determined, the 
mathematical model based on a plurality of par- 20 
tides (13, 14, 15) in a thermal bath (10), in which 
particle(s) (13, 14, 15) at a given location in the 
thermal bath (1 0) correspond to a given base sta- 
tion, and the particle(s) (13, 14, 15) have a plur- 
ality of states, each occupied state correspond- 25 
ing to a radio channel, the occupancy of states 
determined by a thermal distribution, and the 
particles having interaction energies correspond- 
ing to interference between base stations (2, 3, 
4), a desired allocation of radio channels being 30 
determined by setting the initial temperature of 
the thermal bath (10) at infinity and then reduc- 
ing the temperature of the thermal bath in steps 
so that a local minimum energy is attained by a 
stochastic process in which the interaction be- 35 
tween particles (13, 14, 15) corresponds to an 
acceptable level of interference between base 
stations (2, 3, 4) and determining the radio chan- 
nel allocation from the distribution of occupied 
states of the particle(s). 40 

17. A radio system (1) as claimed in claim 16 further 
characterised in that 

- the states of the particles (13, 14, 15) are 

spin states, 45 

- each base station (2, 3, 4) is represented by 
a group of particles (13, 14, 15), 

- the number of particles in a group corre- 
sponding to the number of channels as- 
signed to the base station represented by 50 
that group of particles, and 

- the spin state of a given particle in a group 
represents a particular radio channel. 
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